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Scaling laws of passive tracer dispersion in the turbulent surface layer. 
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Experimental results for passive tracer dispersion in the turbulent surface layer under stable 
conditions are presented. In this case, the dispersion of tracer particles is determined by the interplay 
of three mechanisms: relative dispersion (celebrated Richardson's mechanism), shear dispersion 
(particle separation due to variation of the mean velocity field) and specific surface-layer dispersion 
(induced by the gradient of the energy dissipation rate in the turbulent surface layer). The latter 
mechanism results in the rather slow (ballistic) law for the mean squared particle separation. Based 
on a simplified Langevin equation for particle separation we found that the ballistic regime always 
dominates at large times. This conclusion is supported by our extensive atmospheric observations. 
Exit-time statistics are derived from the experimental dataset and show a reasonable match with the 
simple dimensional asymptotes for different mechanisms of tracer dispersion, as well as predictions 
of the multifractal model and experimental data from other sources. 
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I. INTRODUCTION 

The phenomenon of scalar turbulence or random ad- 
vection of tracer particles by random turbulent flow has 
become a topic of significant attention during the last few 
years Q, 0, 0], [d], [5], [6], Y]. This is not only due to 
its significance for understanding the transport processes 
in global geophysical systems (e.g. the atmosphere and 
oceans) and as a theoretical framework for the design of 
some technological applications (mixers, chemical reac- 
tors, combustion chambers), but is perhaps even more 
important as an instructive example of modern methods 
of theoretical physics applied to a highly non-equilibrium 
dynamic system in order to deduce new phenomenology 
and a wealth of analytical results that can be validated 
numerically and experimentally. Known examples are: 
analytical scaling for white-noise scalar turbulence (solu- 
tions of the Kraichnan model [l|, @], @j); application 
of conformal invariance to the two-dimensional tracer 
flow [J|; re- normalization group formalism; multi-fractal 
structure of tracer statistics 1] and others (see review Q 
and references therein). Remarkably, the recent appli- 
cations of the theoretical framework for this phenomena 
provide rigorous ways to overcome the initially restrictive 
assumptions of the underlying model of locally isotropic 
turbulence by incorporating effects of anisotropy, flow 
boundaries, mean velocity shear, buoyancy, etc. and en- 
able the analytical calculations of associated corrections 

In this paper we report our experimental results on 
the dispersion of passive tracers in the turbulent surface 
layer. The wall-bounded turbulent flow provides flexible 
settings to study the effects of mean velocity shear (i.e. 
non-uniform wind) and system boundaries (the underly- 
ing surface) on the passive scalar dispersion. 



The celebrated Richardson law [111 established the 
growth of the mean inter-particle distance: 



{R^it)) - \e 



(1) 



where n = 3, A is a scale- independent dimensional pa- 
rameter, and became a signature of turbulent dispersion 
(see review [l^l)- It was later recognized by Oboukhov 
and Corsin that this law is a direct consequence of Kol- 
mogorov scaling in turbulence [13|, [3|. Indeed, from a 
power-law assumption for the velocity differences 5v ^ 
dR/dt ex i?'* and (HI) it is easily to derive R ex t^^^ or 
h — 1 — 2/n. The latter expression leads to the Richard- 
son law with {R^{t)) ex t^ for the Kolmogorov scaling 
h = 1/3. This is a reflection of an intimate and well- 
known connection between the power exponent of the 
particle separation law n in ([T]) and the passive tracer 
statistics J3J. From a mathematical point of view, this 
connection can be translated into a scaling law of the 
correlation functions of the tracer concentration that in- 
cludes the parameter n [lj|, 0: 

S2{R) = 2[^2(0) - F2{R)] ex i?2/", n ^ 1/(1 - h),{2) 



where 



F2(i?) = (C(x + R,i)C(x,i)) 



(3) 
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is the pair-correlation function. 

Different values of n (the scaling exponent of the 
mean squared displacement) have been derived theoret- 
ically and experimentally for different kinds of turbu- 
lent flow (see [I4I, and refs therein). Contrarily, a 
particular value of n can be associated with a particu- 
lar energy injection mechanism of turbulence and can be 
used for characterization of its dispersive properties. As 
was mentioned above, for the tracer dispersion by Kol- 
mogorov (locally isotropic) turbulence, n — i (Richard- 
son law). For Bolgiano-Obukhov (buoyancy dominated) 
turbulence, n = 5 [l3|- For turbulence with mean ve- 
locity shear it was recently deduced Q that n = 6 (for 



particle separation along the mean velocity) and n = 4 
(for separation in the transverse direction). It is worth 
noting that in all latter cases the separation of particles 
is always faster than in the standard Richardson regime 
(i.e. without velocity shear). For surface layer turbulence 
(i.e. turbulent motion near a flow boundary), whose scal- 
ing properties are significantly different from the isotropic 
Kolmogorov model, the value of n can eventually devi- 
ate from Richardson's prediction. Particularly, the self- 
similarity arguments applying to this case lead to much 
lower values of the scaling exponent of the mean squared 
displacement (n = 2 or so-called "ballistic" regime, since 
Rent) [1J|. Finally, for the "combined" case (turbulent 
flow near a boundary and a mean velocity shear) Ref Q 
argues that there is no universal value for the parameter 
h in ^ (and hence neither for exponent n). 

From simple dimensional analysis it is obvious that the 
dimensional parameter A in ([1]) should also be different 
for the different values of n, i.e. for the particular type of 
turbulent flow. For instance, in the case of Kolmogorov 
turbulence A = e (e is the kinetic energy flux); but for 
surface layer turbulence, the friction velocity w* (and not 
the dissipation rate e) becomes the only similarity param- 



eter hj, A 



(therefore the expression Xt = (v^,t) 



provides the dimension of i?^). 

The same dimensional analysis also leads to the im- 
portant conclusion that the ballistic regime (and not the 
shear effect) is the universal long-time asymptote of the 
dispersion process of the tracer particle in the turbulent 
surface layer (see [iJI and refs therein) . This can be eas- 
ily seen based on the following arguments. Let us assume 
that in addition to w* , the problem is characterized by a 
length scale I (i.e. initial separation of the particles or the 
initial distance to the ground) . Then the dimensional ar- 
guments translate ([1]) into 



{R^{t))^\t\ \^v:f 



(4) 



A naive conjecture at this point would be that a disper- 
sion mechanism with the highest value of n (e.g. shear) 
will dominate over large time periods (i.e. at i > l/v^). 
However, more careful analysis casts doubt on this con- 
clusion. Indeed, for t » l/v^, particles should "forget" 
about their initial positions and the parameter I should 
be dropped from the expression for the particle separa- 
tion dU all together. We can see that the only possibility 
for this is the case when n = 2: conversion to the ballistic 
regime. 

This variety of possible scenarios of tracer dispersion in 
the surface layer turbulence and some ambiguity of avail- 
able analytical predictions (at least when being straight- 
forwardly applied) motivated our experimental study of 
this phenomenon. 



II. EXPERIMENTAL PROCEDURE 

Scrutinizing measurements at a given location within 
the Taylor hypothesis of frozen turbulence R = Ut, where 



U = const; it follows from ([T]) that: 

S2{R) = 2[F2(0) - F2{R)] oc i?2/" oc (C/i)2/", (5) 

so single point time-measurements of concentration also 
show power-law behavior with the same characteristic 
exponents as the spatial measurements. 

We estimated the value of parameter n in ([5]) based 
on extended atmospheric observations of the tracer con- 
centrations. A continuous (and relatively stable) influx 
of particles into the turbulent flow was supplied by the 
anthropogenic activity in the surrounding urban areas 
(about 20 km in size). This "highly distributed" source of 
particles maintained the quasi-equilibrium (and on aver- 
age, well mixed) tracer distribution which was validated 
by our long-term observations. The observation tower 
was at a height oi H = 12m and the tracer concentration 
was measured by means of scattered light intensity off 
incoming particles. Continuous observation was under- 
taken over a period of 21 days with an air sample taken 
every 5 seconds (30 liters per sample) . The instrumenta- 
tion provided total particle count as well as their size dis- 
tribution in the range of 1 - 10 iim. Local meteorological 
observations (three conventional meteorological stations 
in the vicinity) and global observations (with resolution 
about 1 km) were available to draw the conclusions about 
local meteorological conditions of the surface layer (sta- 
bility, proflles of wind speed and temperature). Wind 
speed and temperature were also measured directly at 
the sampling point {H = 12m). 

As usual in meteorological studies, we anticipated that 
the night time observations would correspond to the case 
of the boundary layer under stable or neutral conditions. 
For convective turbulence to be generated, an intensive 
heat flux from the underlying surface has to occur and 
this, of course, should include the signiflcant effect of 
solar radiation. Indeed, in our experimental results the 
latter condition is bound to the daytime observations. 



III. DATA ANALYSIS 

Undertaking correlation analysis of the time series of 
the concentration C{t), of single point measurements, we 
found that the minimum time span that provides reliable 
statistics of the process should be more than three hours. 
As the flrst step, we plotted ^2 {R) as a function of i? = 
Ut on a log-log scale and evaluated the scaling exponent 
n from the expression 5*2 (i?) oc (Ut)^^"^ predicted by ^ 
and ([5]). Some examples of these plots are depicted in 
Figffl 

In this paper we report on the analysis of the data 
that corresponds to the stable conditions of the turbu- 
lent surface layer (results on convective tracer dispersion 
will be published elsewhere). As was mentioned above, 
the turbulent motion in this case is determined only by 
a one dimensional parameter - the friction velocity u*, 
which provides the scale for velocity fluctuations. This 
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TABLE I. Apropos mean wind speeds Uh for Fig. [T] 
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FIG. 1. Structure function of the tracer concentration 
52 (-R) oc R" oc UHt" on a log-log scale. Inset shows the 
estimation of how the scaling exponent of the mean inter par- 
ticle displacement {R^{t)) oc t" converges over time. Each 
data set corresponds to a maximum of 4 hours of observa- 
tions and wind velocities outlined in Table [T] -f, Sample 1; o, 
Sample 2; *, Sample 3; ■, Sample 4; o, Sample 5; D, Sample 
6; 0, Sample 7; A, Sample 8 and V, Sample 9. a = 2/3 is the 
Richardson regime, a « 2/5 is the velocity shear mechanism 
and a = 1 is the ballistic regime. Global convergence to the 
ballistic regime {a = 1) is clearly visible. 



self-similarity property of the surface layer turbulence is 
expressed by the well-known relations [l3| 



dU{z) 
dz 



Kz' 



Kz' 



(6) 



where z is the vertical distance from the underlying sur- 
face and e is the dissipation rate, k = 0.4. 

The values of wind speed relevant to the data series 
depicted in Fig. [1] are listed in Table HI By using the 
logarithmic profile U{z) = {v^,/K)ln{z/zo) followed from 
^ we can estimate w, = cUh (c = K/ln(iJ/zo)) and 
then e — vI/kzq- Assuming a typical value of roughness 
height zq sa 0.1 m leads to to the estimates c sa 0.01 and 
e w 0.004 kg ni^/s^ for Uh = 5.52 m/s. The parameter 
Th ~ H/U{H) ~ H/cv* corresponds to the time scale 
for a particle to reach the underlying surface and was of 
order 10^ seconds. 

The low relative fluctuations of the wind speed was 
one of the reasons that these data series were selected for 
analysis. It was also a justification for the application of 
the Taylor hypothesis to these datasets. It is worth not- 
ing that a particular value of the mean wind velocity has 
no relevance to the main reported result of our paper; it 
will only change the point at which a data series will ap- 
proach the global asymptote in Fig. [TJ It will not change 
either the existence of, nor the slope of this asymptote. 

The dynamics of particle separation in the turbulence 
surface layer can be described by modifying the Langevin 
equation for particle separation [15;]. By including the 



Sample 



Uh (m/s) 



1. 


4.15 ±0.12 


2. 


2.94 ±0.11 


3. 


2.25 ±0.30 


4. 


3.85 ±0.12 


5. 


5.52 ±0.16 


6. 


7.42 ± 0.05 


7. 


4.27 ±0.21 


8. 


2.49 ± 0.06 


9. 


3.20 ±0.12 



shear term proposed in [&\ we can arrive at the following 
system: 

dR^ = G{R^)dt + D{R)dt ± ^/2K dW, (7) 

dRi_ = D{R)dt ± ^/2K dW, (8) 

where G(Rz) = {dU/dz)Rz is the velocity shear term, 
D{R) = [dK/dR ± 2K/R) is the drift term, K{R) = 
/3e'^^^R*^^ is the diffusivity, /3 = const, e(R) is the local 
dissipation flux (which describes the local activity of the 
turbulent flow) ^ , dW is an isotropic Brownian motion, 
Rz is the vertical component of separation, R^ is its hori- 
zontal component and R^ ~ [Rz , Ry] are the components 
perpendicular to the mean shear; so R^ = R^ + R^. 

For t < Th, U'{z) - v^/kH, e(R) w e{H) = const, 
so G <$:. D. By placing Rz ^ Rx ^ R± ^ R, we recover 
the standard Richardson regime R^ oc t^ [6]. For the 
intermediate times {t w Th) the velocity shear mecha- 
nism dominates (G ^ D) and following arguments of Q 
we arrive at faster separation of particles {R^ ex. i" with 
4 < n < 6). Finally, at the limiting case t > Th the sepa- 
ration becomes slower (ballistic with R? oct^). The latter 
can be deduced from the system ([7]). At f:^ Th we can 
use estimates: U'[z) ~ v^,/kRz, e(R) = v^/kRz, which 
leads to a simplified system for the long-time asymptotes 



/y^R dW, dRi_ ex 



(9) 



since Rz ~ R± ^ R. 

We can see from the diffusive equation that the asso- 
ciated pdf p{R, t) can be derived 



dp 



d_ 
dR 



1V*TT^ I R 



dp 
OR 



(10) 



where 7 is the constant of order of unity. This equation 
has a straightforward solution: 



piR,t) 



1 



27u*t 



exp 



R 

jv^t 



(11) 



It is evident from here that the main asymptote of 
the model as t ^- cx) is the ballistic regime {R{t)) = 
J^ Rp{R,t)dR^vJ. 



The existence of this global ballistic asymptote often 
suppresses the effect of the shear dispersion and clearly 
emerges from our experimental data (see Fig 1). We ob- 
serve that at the short-time intervals, close to i = 0, the 
parameter n « 3 (Richardson regime), then it reaches 
maximum: n « 5 (shear dispersion) and finally n gradu- 
ally decays to the the ballistic value: n k, 2. It is worth 
noting that the long-time ballistic limit n w 2 is in agree- 
ment with recent experimental data on atmospheric dis- 
persion 2 < n < 3, see [ij], [12] and references therein. 

The non-trivial scaling properties of tracer dispersion 
in the turbulent surface layer can be also demonstrated 
by means of exit-time statistics for the concentration time 
series [8| ; which are statistics of time intervals in which 
a measured value of concentration exits through a set 
of thresholds. By scanning the time series for a given 
threshold 5C , one can recover a set of times Ti(5C)^ where 
the measured concentration reaches this threshold. This 
set can then be used to calculate the Inverse Structure 
Function (for details see Q, |10|): 



^,{5C) ^ (r^iSC)). 



(12) 



A comprehensive analysis of the properties of the In- 
verse Structure Function can be fulfilled by applying the 
well-known mutifractal approach |a| , [lal ■ This results in 
the following scaling [8[, |16l |: 



E,(<5C)(x(<5C)^(«), x{q)=mm[{q 

h 



3-D{h))/h] 

(13) 
where h is the index of singularity from the range 
[hmin, hmax] such that t ^ {6C)^ and D{h) is the fractal 
dimension of the set with a singularity h. It is worth not- 
ing that particular values h = 1/3, h = 1/5 and h = 1/2 
correspond to the Richardson law, the shear dispersion 
and the ballistic scaling discussed above. 

According to [I6j, the function {q+3 — D{h))/h reaches 
its minimum at the upper boundary of the singularities' 
range, so we can set h — hmax and write 



X{q) = (g + 3 - D{h,nax))/K 



(14) 



This leads to the important conclusion that xW) is 
a linear function of q (i.e. no intermittency correction), 
which has been verified by numerical simulations and by 
some limited experimental data ^16:i|. This conclusion was 
also validated against our experimental dataset and the 
results are shown in Fig [2] We observe that indeed xil) 
closely follows the predicted linear trend for g > 1 and 
provides a reasonable match with the experimental data 
available in the literature [17|. A change in slope near 
the value q = 1 can be attributed to a contribution of 
the slow (difFerentiable) components of turbulent motion 
[l(j |. A model for this effect will be discussed elsewhere. 

We also present a plot of the mean exit-time Ei((5C) — 
{t{5C)) as a function of the concentration threshold 5C 
in Fig 131 Co and its associated Iq are located at the con- 
centration minima for each sample, and are used to col- 
lapse data only. Three asymptotes corresponding to the 
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FIG. 2. Mean Inverse Structure Function for 12 experimen- 
tal data sets (*). Error bars correspond to ± mean standard 
deviation. The dashed line represents a linear best fit over 
1 < (? < 5 predicted by (|14|) [la]. Experimental values for 
temperature fluctuations (V) of [I3 are presented for refer- 



regimes discussed above are also depicted for compari- 
son. These asymptotes can be easily established from 
([5]) based on simple dimensional arguments. Indeed, 
for the short times t <^ Th we can assume Richard- 
son (or Corsin-Obukhov) scaling SC ex t^^^ and derive 
(t((5C)) ex {6C)^. Similarly, for the shear dispersion we 
can write {t{6C)) oc (^C)", n « 5 Q; while for the ballis- 
tic regime at the longer times {t ^ Th) we arrive at the 
scaling {t{6C)) oc (SC)^. A subsequent change of regime 
as discussed above, corresponding to the curve {t{SC)) 
deviating from one asymptote to another, is clearly vis- 
ible. From this plot it is also evident why in the case 
of wall-bounded turbulent dispersion it is not possible to 
assign any universal value for parameter n in the scal- 
ing law {t{SC)) ex (i5C)" (similarly in ^), which is in 
agreement with the conclusions of 0- 



IV. CONCLUSIONS 

We presented experimental results for passive tracer 
dispersion in the turbulent surface layer under stable con- 
ditions. In this case the dispersion of tracer particles was 
affected by the mean velocity gradient and flow bound- 
aries. We found that our observations can be intrinsically 
explained with the three-stage model of tracer disper- 
sion. During the first stage of separation (T ^ H/v^,, 
where H is the observation height and u* is the friction 
velocity) tracer particles obey the standard Richardson 
model. During the second stage (T « H/v^,) the shear 
mechanism of dispersion dominates [6|. Finally, when 
(T ^ H/v^,) the shear mechanism is followed by a tran- 
sition to the ballistic regime of dispersion induced by the 
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FIG. 3. Exit-time for tiie first moment of the concentration 
time series. Different regimes of dispersion correspond to the 
different power-law asymptotes (see text). Scaling factors Co, 
to were used for collapsing data thus a better visual appear- 
ance (see text). 



specific scaling properties of the turbulence in the surface 
layer. This scenario of inter-particle distance seems to be 
in agreement with atmospheric observations [iJI, uM ^^ 
well as experimental results on tracer dispersion by tur- 
bulent surface flow in water channels [l^ . We found that 
the ballistic regime always suppressed the velocity- shear 
effect predicted in [6] at the later stage of dispersion, 
resulting in much slower rate of particle separation (i.e. 
lower value of parameter n in ([1])). Exit-time statistics 
were derived from the experimental dataset and showed 
a reasonable match with the simple dimensional asymp- 
totes, as well as predictions of the multifractal model and 
experimental data from other sources. 
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